Anderson-Mott Transition in a Magnetic Field: Corrections to Scaling 
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It is shown that the Anderson-Mott metal-insulator transition of paramagnetic, interacting dis- 
ordered electrons in an external magnetic field is in the same universality class as the transition 
from a ferromagnetic metal to a ferromagnetic insulator discussed recently. As a consequence, large 
corrections to scaling exist in the magnetic-field universality class, which have been neglected in 
previous theoretical descriptions. The nature and consequences of these corrections to scaling are 
discussed. 



PACS numbers:71.30.+h; 64.60.Fr; 64.60.Ht 
I. INTRODUCTION 

The metal-insulator transition of interacting, dis- 
ordered electrons at zero temperature (T = 0), or 
Anderson-Mott transition (AMT), has been,studied the- 
oretically near two-dimeasions (d = 2+e) Jad as well as in 
high dimensions (d > 6p and d < 6,ET). In the latter case 
one finds the complete universality that is characteristic 
of Landau theories, while in d = 2+e there is a large num- 
ber of different universality classes. This is because the 
AMT in low dimensions is driven by soft modes, and the 
soft-mode spectrum of electrons at T = is determined 
by the symmetries of the underlying field theory. Accord- 
ingly, systems in external magnetic fields, in the presence 
of magnetic impurities, in the presence of spin-orbit scat- 
tering, etc., all constitute different universality classes. A 
priori it is not clear whether the high-dimensional or the 
low-dimensional theories provide a better description of 
the situation in d = 3. Experimentally, while different 
universality classes seem to exist, the situation is also 
far from clear .a A serious problem is posed by the fact 
that the minimum distance from the critical point, and 
the ease and accuracy with which the distance from the 
critical point can be controlled, are much inferior to what 
can be routinely achieved near thermal phase transitions. 
Consequently, one should expect corrections to scaling to 
play an important role in the experimentally accessible 
region of parameter space, and their incorporation into 
any interpretation of data to be essential. This expec- 
tation notwithstanding, however, corrections to scaling 
near an AMT have been largely ignored by both exper- 
imentalists and theorists, with the exception of a sug- 
gestion that strong (logarithmic) corrections to scaling 
may be responsible for th&UJ-understood observations in 
phosphorus-doped siliconE'El 

Recently, another example of an AMT has been con- 
sidered near two-dimensions, namely, the quantum phase 



transition-from a ferromagnetic metal to a ferromagnetic 
insulator .0 Since both an intrinsic magnetization and an 
external magnetic field break the rotational symmetry in 
spin space and give a mass to the transverse spin-triplet 
particle-hole excitations, one expects this transition to be 
related to the AMT in an external field. An important 
question in this context relates to the Goldstone modes 
or spin waves that are present in a ferromagnet, but not 
in a system in an external field. Since the Goldstone 
modes contribute to the soft-mode spectrum, one might 
expect them to influence the critical behavior. However, 
it has been shown in Ref. |^ that this is not the case. This 
makes it likely that the two transitions indeed belong to 
the same universality class. The effective model for the 
ferromagnetic MIT that was derived and studied in Ref. 
|§|, however, is different from the oneijor the magnetic-field 
universality class considered beforep and the former con- 
tains strong corrections to scaling that were absent in the 
latter. This suggests that the existing description of the 
magnetic-field universality classo is incomplete. 

In this note we show that the two transitions belong 
indeed to the same universality class, as has been sug- 
gested in Ref. ||, and that the previously studied model 
for the magnetic-field AMT missed terms that contribute 
important corrections to scaling. As a result of these cor- 
rections to scaling, a reliable experimental determination 
of the critical exponents is not possible from the currently 
available experimental data. 



II. EFFECTIVE ACTION 

We consider the usual mictoscopic model for an AMT 
in an external magnetic field. cH3 That is, we assume that 
the orbital effects of the magnetic field serve only to sup- 
press thfij soft modes in the particle-particle or Cooper 
channel,B and consider only the Zeeman term explicitly. 
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The latter 's contribution to the action is 

Sb = 1 91, Mb B J dx T £ £ ( x ) ( x ) 

n a 

-^(xX^x)] . (2.1) 

Here t/L is the g-factor, /^b is the Bohr magneton, B is 
the magnetic field, and T is the temperature. ip],i and 
V'nt,!. are fermionic fields for up- and down-spin electrons, 
with n the index of a fermionic Matsubara frequency, 
uj n = 2nT{n + 1/2), and a a replica index to deal with 
the quenched disorder. 

Within the formalism of Ref. |8|, bilinear products of 
fcrmion fields are expressed in terms of a classical matrix 
field Q, whose elements are related to the fermion fields 
by means of the isomorphism 
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ipuihi 
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(2.2) 



Here all fields are understood to be taken at position x, 
and 1 = (rii, cti), etc. In terms of these matrix fields, the 
Zeeman term can be writtenllJ 



Ab = —ih J dx tr [(t 3 <g S3) Q(x)] 



(2.3) 



where h = gi, /^b B, tr is a trace over all discrete indices 
of the matrix Q, and T3 = —S3 = — io~ Zl with o~ z a Pauli 
matrix. More generally, it is useful to expand Q into a 
spin-quaternion basis, 

QM=J2 £ r<9l 2 ( X ) ( T r®Si) , (2.4) 
r=0,3 i=0,3 

with To = so = 1 the unit 2x2 matrix. The partition 
function can be written in terms of a functional integra- 
tion over the field Q and an auxiliary field A,El 



Z= D[Q] D[A] e A[Q - K] , 



with an action 



A[Q,A] = Ab=o[Q,M + A b [Q} 



(2.5a) 



(2.5b) 



A 



B=0i !S 



The action in the absence of a magnetic field, 
the same as in Ref. ||. 

The further formal development proceeds in exact 
analogy to Ref. ||, and we will therefore be very brief. 
The action allows for a saddle-point solution that can be 
written in terms of two Green functions, Q and J 7 , which 
are related to the saddle-point values of the matrix ele- 
ments qQ and 3Q, respectively. They obey the equations 



(*Wn - & - En) £«(k) + A„ T n (k) = 1 , (2.6a) 
(*w„ - a - En) T n {k) + A„ g n (k) = . (2.6b) 

Here £k = k 2 /2m — ep with m the electron mass and 
cf the Fermi energy. The two self-energies E and A are 
given by 



k 

irNpT V ^—^ 



with 



A = ^ + 2r(*)T^l^j-„(k) 



(2.7a) 
(2.7b) 

(2.7c) 



Here is the spin-triplet interaction amplitude defined 
in Ref. ||, r is the transport relaxation time due to the 
quenched disorder, and N-p is the density of state s at the 
Fermi level. The appearance of h in Eq. ( 2.7c ) is the 
only difference between the saddle-point equations here 
and those in Ref. ^, and it results from the Zeeman term 
in the action. For later reference we note that the field 
configuration obtained by putting A = is not a saddle 
point, in constrast to the situation in the absence of an 
external field. 

The discussion of the saddle-point solution proceeds as 
in Ref. |[ In particular, A obeys the equation 



2r<*>A 



V (iuj n - £ k + isgncj„/2r) 2 - A 2 

(2.8a) 

For h — there is a nonzero physical solution only if the 
Stoner criterion NpT^ > 1 is fulfilled. H In paramagnetic 
systems, NpT^ is smaller than the critical value, and 
A = in the absence of a magnetic field. However, for 
nonzero h we have 



A = h 



i + iv F r 



o(b 3 ) 



(2.8b) 



Since A is proportional to the magnetizationja this re- 
flects the fact that a magnetic field polarizes the spins. 

We see that, at the saddle-point level, the only dif- 
ference between the current situation and the one con- 
sidered in Ref. || is the origin of the nonzero value of 
A. More generally, as far as an effective theory that 
concentrates on soft-mode effects is concerned, the only 
difference lies in the existence of Goldstone modes in a 
magnetic system. Since the latter turn out to be irrel- 
evant for describing the metal-insulator transition,El the 
procedure of expanding about the saddle-point, separat- 
ing soft and massive modes, and deriving an effective 
theory for the metal-insulator transition is therefore the 
same in the magnetic field and the ferromagnetic cases, 
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respectively, and we can simply take over the result for 
the latter. The final effective action thus readsu 



^ = — y dx tr [VQ(x)] 2 + 2H J dx tr [SI Q(x)} 
dxtr ((t 3 ®s 3 )[VQ(x)] 2 ) 

(r 3 ® s 3 ) Q(x)l + Ant [Q] ■ (2.9a) 



2G 3 
+2i? 3 / dx tr 



Here 



^12 = (to <8> so) ^12 27rTn , 



(2.9b) 



is a bosonic Matsubara frequency matrix. The bare val- 
ues of the coupling constants G, G 3 , H, and H 3 are 



1/G = 




+ ^0 ) . 


(2.9c) 


1/G a = 




-°"o ) > 


(2.9d) 


i? = 




^ F ) > 


(2.9e) 


#3 = 




-JV?) > 


(2.9f) 



where cr^ and Np are the Boltzmann conductivity and 
the bare density of states at the Fermi level, respectively, 
of systems whose Fermi energy has been shifted by ±A 
from its level in the absence of a magnetic field. No- 
tice that 1/G 3 and H$ vanish for a vanishing magnetic 
field, and for s mall fields are linear in the field, see Eqs. 
( 2.9d , 2.9f , [2.8b| ) . The interaction piece of the action con- 
sists of three terms, 



A int [Q]=4%[Q]+A£ t [Q] + <l[Q] , (2-9g) 



,(t) 



(3) r 



-nT 



K s I dx 6 

Q1Q2 ^ajQrg Oi_2,4 — 3 

•* 1234 

X ^2(-Y tT [(Tr®*0)Ql2(x)] 
r 

xtr [(r r ®s )Q 34 (x)] , (2.9h) 



-nT 



K t / dx 6 

1234 



x^(-) r tr [(T r ®s 3 )Q 12 (x)] 

r 

xtr [(T r (g)s 3 )Q 34 (x)] , 



(2.9i) 



^ 1234 

EE^i'^uW^W , (2-9j) 



rs ij 



where 



in 



7 tr (^Trrl) tr( 



S 3 SjS 



(2.9k) 



and Kt = 2nT^ and K s are the usual spin-triplet and 
spin-singlet interaction amplitudes Ahat are also present 
in the absence of a magnetic fields K3, like I/G3 and 
i/ 3 , is magnetic field dependent, and vanishes linearly 
with the field for small fields. Since it is absent from the 
bare action, but generated by the renormalization group 
at one-loop order ,E3 it is also proportional to the disorder. 
Finally, 



Q12 = Q12 - 812 (r ® so) w ni , 
and Q is subject to the constraints 



Q 2 (x) 



1 



Q f = Q , trQ(x)=0 



(2.91) 
(2.9m) 



III. METAL-INSULATOR TRANSITION 

We recognize Eqs. ( [2.9] ) as the, generalized nonlinear 
<7 model proposed by Finkel'steinEI for the magnetic field 
universality class, except that the latter was missing the 
terms proportional to 1/G 3 , -ff 3 , and K3. As was shown 
in Ref. ||, these terms do not change the asymptotic crit- 
ical behavior. Choosing the correlation length exponent 
v, the critical exponent for the density of states f3, and 
the dynamical critical exponent z as the three indepen- 
dent exponents, we tlj^-Jiave, to lowest order in an ex- 
pansion in e = d — 2,1 



v 

z 



l/e + 0(l) , 
l/2e(l - In 2) 
d , 



(3.1a) 
(3.1b) 
(3.1c) 



and the critical exponent for the conductivity, s = v(d- 
2), is 



s = 1 + 0(e) 



(3.1d) 



While the additional coupling constants 1/G 3 , ff 3 , and 
A" 3 do not contribute to the leading critical behavior, 
they lead to corrections to scaling. As was shown in Ref. 

the least irrelevant operator related to these coupling 
constants has a scale dimension 



A, 



3 In 2 - 2 
1 - In 2 



e + 0(e 2 ) 



(3.2) 



Consequently, in the vicinity of the metal-insulator tran- 
sition any observable u obeys the homogeneity law 



w{t, T,u,...) = b-*»u(ib l/u , Tb z , ub x \ . . .) 



(3.3) 



where t is the dimensionless distance from the critical 
point. Apart from (or instead of) the temperature T, 
other generalized external fields may appear as argu- 
ments of to, e.g., a frequency, or an electric field, depend- 
ing on the nature of the observable under consideration. 
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U denotes the least irrelevant variable (which is a linear 
combination of I/G3, and -K3), is the scale dimen- 
sion of w, b is an arbitrary scale factor, and the ellipses 
denote the dependence of to on variables that are more ir- 
relevant than u. These corrections to scaling were absent 
from the original model of Ref. |[ 



IV. DISCUSSION 

To put our results in context, we give a brief discussion 
of corre ction s to scaling in general, the additional terms 
in Eq. ( [2.9a| ) in particular, and why they were missed in 
previous treatments. Asymptotically close to a cr itica l 
point, observables obey the homogeneity law, Eq. (3.3), 
and of the various coupling constants only those with 
positive scale dimensions need to be kept. Apart from 
the dimensionless distance from the critical point, whose 
scale dimension is by definition the inverse correlation 
length exponent \jv > 0, the only coupling constants 
with positive scale dimensions are generalized external 
fields,E3 like the temperature in our example. Techni- 
cally, the homogeneity law appears as the solution of a 
renormalization group equation, which takes the form 
of a Callan-Symanzik equation or a similar partial dif- 
ferential equation. If one is not asymptotically close to 
the transition, corrections to the [asymptotic scaling be- 
havior appear from two sources :E3 (1) The solution of 
the renormalization group equation becomes more com- 
plicated than a generalized homogeneous function, and 
needs to be expanded in a power series. (2) The "irrel- 
evant operators" , i.e. coupling constants whose .scale di- 
mensions are negative, can no longer be ignored. til Keep- 
ing them, and ordering them with respect-to their irrel- 
evancy, leads to the Wegner expansion.EjO In our case, 
the leading term in the Wegner expansion, which we have 
kept in Eq. (4.1), is more important than the corrections 
to scaling of type (1). In fact, our value — ^A 3 w 0.26 is 
rather small for a leading Wegner exponent, which typi- 
cally is on the order of 0.5. Since these large corrections 
to scaling result fro m the coupling constants I/G3, H3, 
and K3 in Eqs. (2.9), it is important to keep these terms. 

This brings us to our next discussion topic, viz. the 
origin of these additional terms in the action. As we 
have mentioned in Sec. || above, a field configuration 
that has A = 0, and thus describes non-spin polarized 
electrons, is not a saddle-point solution of the field the- 
ory, Eqs. ( |2.5D , underlying the effective action. Refer- 
ence |^ used the implicit assumption that the only effect 
of the magnetic field, as far as the metal-insulator tran- 
sition is concerned, is to give a mass to certain modes 
that are massless in the absence of an external field E If 
this were true, then one could use a field configuration 
with A = as the starting point for deriving an effec- 
tive action. The result of this procedureo is Eqs. ( |2.9| ) 
with I/G3 = H 3 = K 3 = 0. By the nature of the renor- 
malization group flow equations,! zero bare values for 



these coupling constants imply that they are not gener- 
ated under renormalization cither. The effective theory 
that was originally proposed for this problem is therefore 
incomplete. As we have seen, this omission has an in- 
fluence, not on the asymptotic scaling properties, but on 
the corrections to scaling, within the resulting effective 
field theory. Since the additional terms discussed above 
describe a physical effect, namely the spin polarization 
due to the magnetic field that was neglected in the orig- 
inal treatment, we conclude that the leading corrections 
to scaling discussed in this paper are an immediate, if 
a priori not entirely obvious, consequence of this spin 
polarization. We also note that there is no obvious rea- 
son underlying our result that the additional terms do 
not change the asypmtotic critical behavior; this appears 
to be accidental. In particular, it needs to be stressed 
that this conclusion holds only in an e-expansion near 
two-dimensions, and in d = 3 the new terms might well 
influence the asymptotic scaling behavior. 

We finally discuss the experimental implications of our 
results. For definiteness, let us consider the conductivity 
cr, whose scale dimension is x a ~ sjv. A common way to 
analyze experiments is to extrapolate the data to T = 
and consider a(t,T — 0). From Eq. (|3.3[), we have 



(4.1) 



a(t, T = 0) = cr t s [l + const, x t~ vX *] 



Here o~q is a microscopic conductivity scale on the or- 
der of the Boltzmann conductivity, and the constant is 
nonuniversal. Assuming that the constant is of 0(1), in 
order to achieve, e.g., a 10% accuracy in a determination 
of the value of s from a log- log plot, the i-range needs 
to be restricted to t < (0.1/^|A 3 |) 1/l/|A31 . Extrapolat- 
ing the results of our one- loop approximation to d = 3, 
we have f | A3 1 ~ 0.26, or t < 0.026. If the copstant 
is larger, as is the case at many critical points, Ii3 the 
asymptotic critical region is even smaller. For instance, 
for const. = 1/|^A 3 | ss 4, we need t < 10~ 4 to measure 
s with a 10% accuracy. Such small values of t are not 
currently achievable for Anderson-Mott transitions. The 
technique that yields the best control over t, viz. stress 
tuning, allows to probe the region 10~ 3 < t < 10~ 2 , while 
other methods are restricted to t < 10~ 2 or largerE We 
conclude that it is currently not possible to determine 
the asymptotic critical exponents for the magnetic-field 
universality class of the Anderson-Mott transition with 
any meaningful accuracy. 
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